On smoothness of Black Saturns 

Piotr T. Chrusciel* 
--^l Institute of Physics, University of Vienna 

O Michal Eckstein 

Instytut Fizyki, Uniwersytet Jagiellonski, Krakow 



(N 



Oh 






Sebastian J. Szybka"!" 
Obserwatorium Astronomiczne, Uniwersytet Jagiellonski, Krakow 

July 22, 2010 



Abstract 



, We prove smoothness of the domain of outer communications (d.o.c.) 

j>. of the Black Saturn solutions of Elvang and Figueras. We show that 

00 the metric on the d.o.c. extends smoothly across two disjoint event 

^O horizons with topology ]R x S''^ and E x 5^ x S^. We establish sta- 

ble causality of the d.o.c. when the Koniar angular momentum of the 
spherical component of the horizon vanishes, and present numerical 
C ■ evidence for stable causality in general. 
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1 Introduction 

In [ij , Elvang and Figueras introduced a family of vacuum five-dimensional 



asymptotically fiat metrics, to be found in Appendix [A. 1[ and presented ev- 
idence that these metrics describe two-component black holes, with Killing 
horizon topology M x ((S"^ x S'^)US^)) . In this paper we construct extensions 
of the metrics across Killing horizons, with the Killing horizon becoming an 
event horizon in the extended space-time. Now, it is by no means clear that 
those metrics have no singularities within their domains of outer communi- 
cations (d.o.c), and the main purpose of this work is to establish this for 
non-extreme configurations. Again, it is by no means clear that the d.o.c. 's 
of the solutions are well behaved causally. We prove that those d.o.c. 's are 
stably causal when the parameter C2 vanishes (this condition is equivalent 
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to the vanishing of the Komar angular momentum of the spherical compo- 
nent of the horizon, compare [4, Equation (3.39)]), and present numerical 
evidence suggesting that this is true in general. 

Given the analytical and numerical evidence presented here, it appears 
that the Black Saturn metrics describe indeed well behaved black hole space- 
times within the whole range of parameters given by Elvang and Figueras, 
except possibly for the degenerate cases when some parameters Oj coalesce, 
a study of which is left for future work. In particular we have rigorously 
established that the Black Saturn metrics with C2 = and with distinct a^'s 
have a reasonably well behaved neighbourhood of the d.o.c. Our reticence 
here is related to the fact that we have not proved global hyperbolicity of 
the d.o.c, which is often viewed as a desirable property of the domains of 
outer communications of well behaved black holes. In view of our experience 
with the Emperan-Reall metrics [2] , the proof of global hyperbolicity (likely 
to be true) appears to be a difficult task. 

We use the notation of [4j, and throughout this paper we assume that 
the parameters ai occurring in the metric are pairwise distinct, Oj / aj for 
i +3- 

2 Regularity at z = ai, p = 0, and the choice of c\ 

We consider the metric coefficient gtt on the set {p = 0, z < a\\. A Math- 
EMATICA calculation shows that gu is a rational function with denominator 
given by 

- (2(03 - ai)(a2 - 04) + (as - a\)c\C2)^ (z - ax)(z - a2)[z - 04) , (2.1) 

which clearly vanishes as z approaches a\ from below (we will see in SectionW^ 
that the first multiplicative factor is non-zero with our choices of constants). 
On the other hand, its numerator has the following limit as z — >• ai, 

(02 - ai)^(a3 - ai)(a5 - ax) (2(03 - a{){a^ - ax) - (05 - ai)ci) C2 , (2.2) 

which is non-zero unless C2 vanishes or ci is chosen to make the before-last 
factor vanish: 

ci = ±./ '^"^""^^^"^""^^ /0. (2.3) 

V as - ai 

This coincides with Equation (3.7) of pj. 

By inspection, one finds that the metric is invariant under the transfor- 
mation 

(Cl, C2, il) ^ (-C1, -C2, -V') • 
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Thus, an overall change of sign (ci,C2) i— )■ (— ci,— C2) can be implemented 
by a change of orientation of the angle ^l'. Hence, to understand the global 
structure of the associated space-time, it suffices to consider the case 

ci >0; 

this will be assumed throughout the paper from now on. 



If ( 2.3 ) does not hold, the Lorentzian norm squared gtt = g{dt, dt) of the 
Killing vector dt is unbounded as one approaches ai ; a well known argument 
shows that this leads to a geometric singularity. 



We show in Section 5.8.1 that the choice (2.3) is necessary for regularity 
of the metric regardless of whether or not C2 = 0: without this choice, (7^^ 
would be unbounded near ai, leading to a geometric singularity as before. 



With the choice (2.3) of ci, or with C2 = 0, the point ai := (p = 0, 2 = ai) 
in the quotient of the space-time by the action of the isometry group becomes 
a ghost point, in the sense that it has no natural geometric interpretation, 
such as a fixed point of the action, or the end-point of an event horizon. 
Now, the functions 

Ri := yjp^ + {z- ai)2 

are not differentiable at p = 0, 2; = Oj. So, a generic function of Ri will have 
some derivatives blowing up at /o = 0, z = ai. However, this will not happen 
for functions which are smooth functions of i?f . It came as a major surprise 
to us that the choice of ci above, determined by requiring boundedness of 
gtt on the axis near ai, also leads to smoothness of all metric functions near 
z = ai. It turns out that there is a general mechanism which guarantees 
that; this will be discussed elsewhere [3]. 

To establish that the metric is indeed smooth near the ghost point ai, 
we start with 

gtt = -TT = --F^TT ='■ ^(W,/^A,Ci,C2,/) ) , 
Ilx r Ux 

where A runs from two to five. $ is a rational function of its arguments, 
and hence a rational function of Ri . So gtt will be a smooth function of R^ 
near i?i = if and only if $ is even in Ri : 

$(i?i - {z- ai),fj.A,ci,C2,p'^) = $(--Ri - (z-ai),/iA,ci,C2,p^) , (2.4) 

assuming moreover that the right value of ci has been inserted. (We em- 
phasise that neither FH^ or FHy are even in i?|, so there is a non-trivial 
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factorisation involved jj moreover gu is not even in Ri for arbitrary values of 
the Cj's, as is seen by setting ci = C2 = 0.) Now, there is little hope of check- 
ing this identity by hand after all functions have been expressed in terms 
of p, z, and the aj's, and we have not been able to coerce Mathematica to 
deliver the required result in this way either. Instead, to avoid introducing 
new functions or parameters into <I>, we first note that 

-Ri - {z- ai) = , 



and so (2.4) reads 



$(//!, //A, Cl,C2,p^) = $( ,fiA,Ci,C2,p'^) . 

Ml 

From the explicit form of the functions FH^ and FHy we can write 

where the $j's are polynomials in cf, pi and p^, and G is a polynomial in 
Pi, ci, C2 and p^. One then checks with Mathematica that each of the 
coefficients $j has a multiplicative factor that vanishes after applying the 



identity (5.1) below to replace each occurrence of cj in terms of the pis: 



2 ^ (-/^l + P^){-Pl + P4)P5{f^lP3 + p^){piPi + p^) 
^ PlP^P4{-Pl+ Pb){lJ'll^b + P^) 

It is rather fortunate that each of those coefficients has a vanishing factor, as 
we have not been able to convince Mathematica to carry out a brute-force 
calculation on all coefficients at once. 

An identical analysis applies to gpp = Qzz and u)^/Hy; regularity of 
g^^ immediately follows; there is nothing to do for g^,^. Before doing these 
calculation, care has to be taken to eliminate, with the right signs, all square 
roots of squares that appear in the definition of oj^ . 

3 Asymptotics at infinity: the choice of q and k 

We wish to check that the Black Saturn metric is asymptotically flat. As 
a guiding principle, the Minkowski metric on M^ is written in coordinates 



^We are grateful to H. Elvang and P. Figueras for drawing our attention to the fact 
that this factorisation takes place in the Emparan-Reall limit of the Black Saturn metric. 
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adapted to U(l) x U(l) symmetry as 
V = 



dt^ + dx^ + dy^ + dx^ + djf 



-dt^ + dp^ + p^di? + dp^ + p^dLp^ , 



(3.1) 



with 



(x,y) = p(cos^,sin'(/') , {x-,y) = /5(cosc/?,sin(/?) . 
Introducing p and 6 as polar coordinates in the (p, p) plane, 

(p, p) = r(cos ^, sin 9) , 



the metric (3.1) becomes 

T] = -dt^ + dr^ + r^d6'^ + r^ sin^ 61 dV'^ + r^ cos^ 6* d(^^ 



(3.2) 



Note that 9 G [0, 7r/2] since both p and p are positive in our range of interest. 
As outlined by Elvang and Figueras in W, relating the (p, z, "0, V') coor- 



dinates of the Black Saturn metric to the (r, 0, -0, (p) coordinates of (3.2 ) via 
the formulae 



1 2 1 2 

p = -r sin 2^ , z = -r cos 20 , 9 G 



0, 



vr 



(3.3) 



should lead to a metric which is asymptotically flat. Under (3.3) the metric 



(3.2) becomes 

T] = -dt^ + r''^{dp'^ + dz^) + r^ sin^ 9 dip"^ + r^ cos^ 9 dip^ 



(3.4) 



so that in such coordinates a set of necessary conditions for asymptotic 
flatness reads 



gtt 



-1 , r ^sin ^6'5rtv -^ , 



r'^Qpp = r'^Qzz -> 1 , r ^ gin ^ g^^ -^ 1 , r ^ cos" 



'5w 



(3.5) 
1, (3.6) 



when r tends to infinity. One also needs to check that all metric components 
are suitably behaved when transformed to the coordinates (x, y, x, y) above. 
Finally, each derivative of any metric components should decay one order 
faster than the preceding one. 
We start by noting that 



1 



-r cos 



sin^e) = -(p2 



f) 
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which is a smooth function of (x, y, x, y). On the other hand, 

p = r sin 9 cos = pp 

is not smooth, but its square is. This imphes that all the functions appearing 
in the metric are smooth functions of {x,y,x,y), except perhaps at zeros of 
the functions Ri and of the denominators; the former clearly do not occur at 
sufficiently large distances, while the denominators have no zeros for p > 



by Section 5.3, and at p = away from the points Cj by Sections 5.4 and 

Mm 

To control the asymptotics we note that pi = 0{r'^), but more precise 
control is needed. Setting R^ := p^ + z^ = r^/4, a Taylor expansion within 
the square root gives 



l^i 



\/p2 + (2; - aj)2 - {z- ai) 
[ 2zai - af 



a. 



2 



= (r^ + 2a, + 2^(1 + cos 2^)) sin^ 6 + 0{r-^) 
= {r'^ + 2ai)sin^e + 0{r-^) . 

For z < this can be rewritten as 

Pi = {r'^ + 2ai + 0{r~^))sin^e . (3.7) 

To see that the last equation remains valid for z > we write instead 



A^i 



P' 



sj p^ ^{z- aiY + (z - ai) 
i?2 sin2 2d 



RyJl-^-^^ + {z-ai) 
Rsm'^29 

Rsm'^29 
(l + |)(l-f + 0(i?-2)) 

1 + cos 29 



and we have recovered (3.7) for all z, for r large, uniformly in 9. 
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The above shows that fj,i — fij = 0{1) for large r; in fact, for i j^ j, 
fii - fij = (2(aj - Qj) + 0{r^'^)) sin^ 9 . 



Keeping in mind that 



p + fiifij ~ r sin I 



where we use / ~ g to denote that C ^ < f/g < C for large r, for some 
positive constant C, we are led to the following uniform estimates 



Mo « r^o sin^s i 
Ml f« r24 sin^s e sin^ 26* , Mi ^e!_ 
M2«r28sin24 6lsin2 26l, M2 



r24sin24^sin4 2^ 



M3 ss r^o sin^'^ 



M1M2 

i^«r28sin28| 



,30 ^;„26 , 



, M4 « r^" sin^ 
F ^ r^s s^j^34 Q ^ 

a. = . . o . (1 + 0{r 2)) « r^ -"-2 



4 sin^ I 

^2\2/ 



cos 

2\ ^ ^16^:^8, 



P = (/"3 /W4 + P^)^(/Ui Ai5 + P^)(/W4 ^5 + P^) ~ ^'^^ sin 
This shows that, for large r. 



H:r 



F- 



28^-28 , 



Mo + ci C2 Ms + 44 M4 +0(r^« sin^« 9) 



H,, 



F 






Mo^ + ci C2 M3 + c\cl mA +0{r^'' sin^^ , 

P2 jJ-l 



_30 <,;,i26 ( 



and in fact the ratio tends to 1 at infinity. We conclude that 

gtt + l = Oir-^) , 

uniformly in angles. 

In order to check the derivative estimates required for the usual notion 
of asymptotic flatness, we note the formulae 

Pi = m + l/2( -x'^ -f + x^ + f 



+ jAai — 4aj(x2 + y"^ — x^ — y^) + (x^ + y^ + x^ + y^) 
{x^+f){x^ + f) . 
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Since the //j's and p^ are smooth functions at sufficiently large distances, 
it should be clear that every derivative of any metric function decays one 
power of V x^ + y^ + x^ + y^ faster than the immediately preceding one, as 
required. 

The constant q appearing in the metric is determined by requiring that 
fft^ — >• as r tends to infinity. Equivalently, since gu — >• —1, 

a = — lim 



r->oo Hy 



Now, 



UV _ ^ Ci i^i VMqMi - C2 R2 VM0M2 + cf C2 R2 VM1M4 - Ci 4 Ri VM2M4 

Hy FHy^/G^ 

H4 R2VMM + ciC2RiVIhm + 0{r^^) 
''V3VG;(Afo^ + ciC2M3 + c2ciM4f + 0(r28)) ' 



M2 " "■ -^ ■ i z -'111 

where we have not indicated the angular dependence of the subleading terms, 
but it is easy to check that the terms kept dominate likewise near the axes. 
A Mathematica calculation gives 

_ 2C2KI 

2ki — 2kiK2 + C1C2K3 

which can be seen to be consistent with |4|, when the required values of the 
Cd's are inserted. 



In view of (3.6), the constant A; > needs to be chosen so that 



k^ lim r^H^P = 1 . 



One finds 

4Kf(-l + K2) 



k' 



-2k,i{-1 + K2) + CiC2K3)2 



as in |4j. From (3.7) and from what has been said so far one immediately 
finds 

HxCry 



lim r ^ sin '^6 qm = hm ^ . o^tt 

r-s>oo ^^ r^oo r^ SU\ OH, 



y 



r Gy ^3/i5 

hm -— — ^-K— = lim 



r-s>oo 7-2 sin^ 9 r^°° r^ sin^ 6 //4 
1, 
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Gx ,. P^/i4 



as desired. Finally, it is straightforward that 

hm r"^cos"^6'5(p<p = lim -^ — %-^ = hm 

= 1 . 

Further derivative estimates follow as before, and thus we have proved: 

gpLv - V^y = 0{r'^) , di^... di^g^y = 0{r'^~^) . (3i 

4 Conical singularities and the choice of C2 



It is seen in Table 5.1 below that g^^p vanishes for {z < as} U {04 < z < 03}, 



while gpp does not, which implies that the set {z < 05} U {04 < 2; < 03} is 
an axis of rotation for 9^. In such cases the ratio 

hm^ 

determines the periodicity of ip needed to avoid a conical singularity at 
zeros of d^, and thus this ratio should be constant throughout this set. This 
leads to two equations. For {z < oi}, the choice of k already imposed by 
asymptotic flatness leads to 

lim ^p^ = 1 . (4.1) 

P-^0 dipip 

Either by a direct calculation, or invoking analyticity at p = across z = a^, 
one finds that the same limit is obtained for ai < z < a^ with the choices 



of k ad ci determined so far. The requirement that (4.1) holds as well for 
a4 < z < 03, together with the choice of k already made, gives an equation 
that determines C2: 

lim -^p = 2{a2 — ai)(a3 — 04) x 

/ (03 - ai)(a2 -Q4) 

Y (02 - 05)(«3 - a5)(2(a3 - ai)(a2 - 04) + (05 - ai)ciC2)2 
= 1 . 

Therefore, to avoid a conical singularity one has to choose 

(03 - ai)ci5i ± (ai - a2)(a3 - a4)5'2 ,. „^ 

C2 = ^ -, TT 77 ^^ , (4.^j 

(ai - a5){a5 - 02) (05 - a3)cf 
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where 

5i = (a2 - a4)(a2 -05)(a3 - as) , 
S2 = y (as - ai)cf5i . 



Equivalently, 
C2 = v2(a4— a2 



±(ai - Q2)(Q3 -Q4) + \/(ai - a3)(a4 -Q2)(Q2 - a5)(a3 - Q5) 
Y^(ai - a4)(o2 - a4)(ai - a5)(a2 - 05)(a3 - 05) 

as found in [4j. 

The case C2 = 0, which arose in Section [2| is compatible with this equa- 
tion for some ranges of parameters Oj, we return to this question in Sec- 
tion ESU 

It follows from the analysis of Section [3] that the analogous regular-axis 
condition for z > 02, 

is satisfied at sufficiently large distances when k assumes the value deter- 



mined there. One checks by a direct calculation (compare (5.30)) that the 



left-hand side of (4.3) is constant on (02,00), and smoothness of the metric 



across {p = 0, z G (02,00)} ensues. 

5 The analysis 

5.1 The sign of the /ij's 

Straightforward algebra leads to the identity, for i ^ j, 

_ (p-i - fJ'k){p'^ + PifJ'k) /^ -,^ 

O'i ~ Uk — - . W-J^J 

Since all the /ij's are non-negative, vanishing only on a subset of the axis 

^:={p = 0}, 
we conclude that 

the Pi — /ifc's have the same sign as the Oj — a^'s. (5-2) 



Furthermore from (5.1) we find 



ai+2-ai {pi+2 - pi){p + piPi+i) ^ f. .^ „s 

Ki := = ^ > . (5.3) 

02-01 2pipi+2\a2 - ai) 

We infer that the functions M^, v = 0, ... ,4 are non-negative: indeed. 



this follows from the fact that the /Xj^'s are non-negative, together with (5.2) 
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5.2 Positivity of H^ for p > 

We wish to show that Hx is non-negative, vanishing at most on the axis 
{p = 0}; note that by the analysis in Section^ Hx certainly vanishes at 
9 = 0. 

Now, Hx vanishes if and only if its numerator vanishes: 

Mo + cJMi + C1M3C2 + (M2 + clM4)cl = . (5.4) 

This equation may be seen as a quadratic equation for C2; its discriminant 

A = cJMl - 4{Mo + cJMi){M2 + cJM^) 

can be brought, using Mathematica, to the form 

A = -4(/ii -/i2)^M2/^3(At2 -^4)^/^4/U5P^(/Wl/^2 + P^)^(/U2/^3 + P^)^ 

X (/X2Ai5 + P^f (ciPi/X3/i4(/ii - ^5)^ - (^1 - ^3)^/^5(^1/^4 + P^fj 

< , (5.5) 

the last inequality being a consequence of the non- negativity of the /ij's. 
Therefore, if a real root exists away from the axis £/, then A = at the 
root and cf satisfies there 



(/^i -/^3)^/^5(/^l/^4 + P ^^^ 



Ci = 2 ? ^2 • <-^-6) 

//f//3/^4(/^i - i^5r 

On the other hand, the smoothness of the metric at /> = implies (com- 



pare (2.3)) 



2 _ r22KlK2 



,_^^ (5.7) 

where, following [4], L is a scale factor chosen to be L^ = 02 — ai. We rewrite 



(5.7) with the help of (5.3), 



^2 ^ ip3 - fJ-i)iPi - /^l)/^5(/^l/^3 + P^)(/^i/^4 + p"^) ,^ 

^ PlP3P4.{P5 - P-lliPlPb + P^) 



Subtracting (5.6) from (5.8) leads to the equation 



{p-i - P3)Pb{Pl + P ){PlPA + P ) .. 
P'lp-iPiiP'i - At5)^(/^i/^5 + P^) 

i^fllfl-iifll - /i4) + pi(/X4 - P3)p5 + (pi - p3)p^j = . (5.9) 



It follows from (A.15), (5.2), and from non- negativity of pj that each term 



in the last line of (5.9) is strictly negative away from £/. We conclude that 



this equation can only be satisfied for p = 0, hence Hx is non-zero for p 7^ 0. 
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5.3 Regularity for p > 

In this section we wish to prove that the Black Saturn metrics are regular 
away from the axis p = 0. For this it is convenient to review the three- 



soliton construction in [|4!|. The metric (A.l) was obtained by a "three- 
soliton transformation" , a rescaling, and a redefinition of the coordinatesjj 
The following generating matrix 

^Q{\,p,z) = diag<- -,- — TT.-T^ TT f 5.10 

t(/i4-A) (/i2-A)(p5- A) (P5-A)J 

was used, starting with the seed solution 

Go = diag{J.,^,-^|, (5.11) 

1/^4 /U2Ai5 Pb) 

where p.^ = —p^/fi^. The general n-soliton transformation yields a new 
solution G with components 

^ /^ N >^ {GQ)acmf\T-^)kimf {GQ)db ., ^„. 

^ab = \^0)ab - > , ^^; l^-l^J 

(the repeated indices a,b,c,d = 1, . . . ,D — 2 are summed over). The com- 
ponents of the vectors m^ ' are 

mi'^=mi'^[%Hiik,P,z)]^^, (5.13) 

where w-q^ are the "BZ parameters" . The symmetric matrix T is defined as 

r.. = "^"''j;°j°^,"^'' , (5.14) 

P^ + PfcW 



and the inverse T ^ of F appears in (5.12). Here pi stands for pi for those 



i's which correspond to solitons, or pi for the antisolitons, where 

Pi = -\/ p^ + (2 - flj)^ - (^ - «i) • 
The three-soliton transformation is performed in steps: 

• Add an anti-soliton at 2: = ai (pole at \ = p\) with BZ vector ttiq 
(1,0, ci), 



■^It has been mentioned at the end of Sec. 2.2 of [4] that the same solution can also be 
obtained (in a slightly different form) as a result of a (simpler) two soliton transformation. 
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• add a soliton at z = a2 (pole at A = ^2) with BZ vector ttiq = 
(1,0,02), and 

• add an anti-soliton at z = 03 (pole at A = /X3) with BZ vector rrig = 
(1,0,0). 

Recall the ordering ai < 05 < 04 < 03 < 02, and we impose the regularity 



condition (5.7). Using these assumptions, we show that that the procedure 



described above leads to a smooth Lorentzian metric on {p > 0}. 
Firstly, we note that 

• ^j — /ifc / for i 7^ /c and p > 0, 

• /ij - /ifc / for p > 0, 



where the first point follows from (5.1). The second statement is a conse- 

^ + \/p^ + {z - o-kY + Q^ - Q fc, hence 

Cfc)^)^, which 



quence of: ^t — fik = \/ P^ + (z — a 

Pi-/2fc = implies (oj-Ofc)^ = (a/p^ + (^ - ai)^ + \//'^ + ( 

is equivalent to 



(? + v^p2 + (2; - aiY^J p^ + (2 - afc)^ + (^ - «j)(^ - Ofe) = . 

The middle term dominates the absolute value of the last one, which implies 
that the last equality is satisfied if and only if p = and (z — Cj) (z — afc) < 0, 
in particular it cannot hold for p > 0. 

We conclude that -0^ is analytic in p and 2; on {p > |. Su bsequently 
the components of the vectors m!^ are analytic there (see (5.13)) and so is 
the matrix F (see (5.14)). The n-soliton transformation (5.12) contains F~^, 



thus detF appears in denominator in all terms in sum in (5.12) (excluding 
(G'o)afe)- Since the numerator of these terms contains analytic expressions 
and a cofactor of F, then only the vanishing of det F may lead to singularities 
in the metric coefficients gab on {p > O}. We show below that detF does 
not have zeros there provided that the free parameters satisfy the regularity 



conditions (5.7). This will prove that the metric functions gtt-, Qt-ip and g^^ 
are smooth away from {p = 0}. Hence 



Hx 



Hx 



Hy 
Hx 



Ho, 



+ <lf 



GyHx 

H„ 



are smooth for p > 0. This is equivalent to smoothness, away from the axis, 
of the set of functions 



Hy 

Her. 



Hx 



w; 



■^ 



HyHx 
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Since H^ has been shown to have no zeros away from the axis, we also 
conclude that 

Hy 

is smooth away from p = 0. 

The next steps in the construction of the line element (A.l) involve a 
rescaling by p^-^^^ and a change of t, ^ coordinates i — )■ t — q"^, ^ — t- — ^. 
These operations do not affect the regularity of the metric functions. 

Let us now pass to the analysis of det T. The metric functions Qpp = Qz-, 
denoted as e^'^ in |4j, can be calculated using a formula of Pomeransky 10 



H^k'P^e''^ = e''^^^^^, (5.15) 

where [4] 

21.0 _ t.2 /^2 ;^5(P^ + m ^J'2Y{P^ + ^J'l /^4)(P^ + mi At5)(P^ + M2;^3)(/9^ + /^3 ;^4)^(P^ + M4 pts) 
/^l {P^ + /^3M5) (/5^ + /^l Ms) (P^ + M2 M4) (P^ + M2 M5) ni=l (P^ + Pi ) 

(5.16) 



and where r^*^-* corresponds to T with ci = C2 = 0. But from what has been 

said the functions detr*-*^' and P do not have zeros for p > 0. Since we have 

shown that H^ does not have zeros there, the non- vanishing of det F follows. 

We conclude that the metric functions appearing in the Black Saturn 



metric (A.l) are analytic for p > 0. It remains to check that the resulting 
matrix has Lorentzian signature. This is clear at large distances by the 
asymptotic analysis of the metric in Section [Sj so the signature will have the 
right value if and only if the determinant of the metric has no zeros. This 
determinant equals 

det 5^^ = -p^Hlk^P"^ . (5.17) 



and its non- vanishing for p > follows from Section 5.2 



5.4 The "axis" {p = 0} 

The regularity of the metric functions on the axis {p = 0} requires separate 
attention. The behaviour, near that axis, of the functions that determine 
the metric depends strongly on the part of the z axis which is approached. 
For example, the pj's are identically zero for z > aj at p = 0, but are not 
for z < tti. This results in an intricate behaviour of the functions involved. 



as illustrated by Tables 5.1 and 5.2 
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z 


P 


G^ - 9.. - ^:;:: 


/-< _ M3M5 
•^S/- M4 


z<ai 


2\z-a^f{z-a,f(z-ai)(z-a,f 


2-04 2 


2(2-a3)(2-a5) 
2-a4 


Ol < z < 05 


2^[z-a,f{z-a,f{z-a,)[^iE^)p' 


z-ai 2 
2(2-a3)(2-a5)/^ 


2(2-a3)(2-a5) 
2—04 


as < z < 04 


2\z-a,)\z-a,f['^)p^ 


2(2-a4)(2-a5) 
z-az 


z-a-i 2 
2(2-a4)(2-a5)^ 


04 < z < 03 


(Sf)%' 


z-a,, 2 
2(2-a4)(2-a3)^ 


2(2-a4)(2-a3) 
2-05 


03 < z < 02 


p' 


2(2-a3)(2-a5) 


2-04 2 
2(2-a3)(2-a5)'^ 


a2 < z 


P' 


2(2-a3)(2-a5) 
2—04 


2-a4 2 



Table 5.1: Leading order behaviour near p = of P, Gx and G 



y 






Hrr 



9wl9, 



pp 



(2(ai-a3)(a2-a4)+(ai-a5)ciC2)^ ~ 

2^^{ax-ao^'^(a2-aAY{ai-z){az-zY{a4-zY{z-ar,Y 



z < a\ 



2^[ax~a?y(a2-aA) „2 _ 2 

(2(ai-a3)(a2-a4)+(ai-a5)ciC2)^fc^ " " 



n. ^ y ^ n. {a2Cx-a\C2^a^{c'i-cx)Y{z-ax) 

ai < Z < as 2«(ai-a3)(ai-a4)(a2-a4)"(a3-^)^(a4-^)^(a5-z)^^ 



2(a3-ai)(a4-ai)(a4-a2) _ 

(a5-ai)(a2Ci-aiC2+a4(c2-ci))^fc^ ^ P 



as < z < 04 



P 



1 (black ring horizon 1\ 



flA < Z < n-i (ai-a2Y{ai-z){z-a^) 

a4 < 2: < a3 2(ai-a3)(a2-a4)(a2-a5)(a3-a5)(a3-^)^ 



(ai-a3)(a2-a4)(a2-a5Ka3-a5) 2 
(ai — a2)'^(a3— Q4)'^A?~ " 



a3 < z < a2 



P" 



1 (spherical horizon ? 



?^ 



a2 < z 



1 



Table 5.2: Leading order behaviour near p = of H^ and of gip^pj Qpp- The 
value 1 of the coefficient in front of f? is precisely what is needed for absence 
of conical singularities at the axis. We write f ^ p"', for some a S M, if the 
leading order behaviour of /, for small p, is f = Cp"', for some constant C 
depending upon the parameters at hand, the exact form of which was too 
long to be displayed here. The question marks concerning the horizons are 
taken care of in Section I5.5H5.7I 
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5.4.1 g^^ 

A complete description of the behaviour of 5^^ at p = can be found in 



Table [57T| One can further see from Table 5.2 that the Killing vector field 
9(p has a smooth axis of rotation on {p = 0, z < as} U {p = 0, 04 < z < 03}, 
as already discussed in Section |4| 

5.4.2 gtt 

At p = 0,z < ai, the metric function gu is a rational function of z with 
denominator 



a{ai — z){a2 — z)(a4 — z) , (5.18) 



where 



a := (2(ai - a3)(a2 -a4) + ciC2(ai -as))^ 

_ 4(ai - 02)^(01 - a3)(a2 - Q4)(Q3 - «4)^ 
(02 - as) (05 - 03) 

So a is nonzero when all the Oj's are distinct. We have already seen that 



the singularity at 2; = oi is removable; the ones suggested by (5.18) at 02 
and 04 are irrelevant at this stage, since we have assumed z < ai to obtain 
the expression. 

From what has been proved in Section [2j gu extends analytically across 
z = ai, so the last analysis applies on p = 0,ai < z < a^. 

The zeros of the denominator of gu restricted to p = 0, a^ < z < a4 turn 
out not to be obvious. It should be clear from the form of gu that those 
arise from the zeros of the numerator of Hx- This numerator turns out to 
be a complicated polynomial in the Oj's, z, and the Cj's, quadratic in C2Jj As 
in Section [2j we calculate the discriminant of this polynomial, which reads 

8(ai - 02)^(01 - ai)^{ai - z)'^{a2 - 04)^(02 - as)^(a3 - z){a4, - z){a5 - z) , 

and which is negative because of the last factor. We conclude that gu does 
not have poles in (as, 04). 

The apparent pole at z = as above is removable: Indeed one can compute 
the limit z ^ a^ using the formula for gu at p = 0, z G (ai, as). After ci is 



^The reader is warned that the numerators listed below depend upon whether or not 
the constants Ca and k have been replaced by their values in terms of the ai's. 
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substituted, one obtains a rational expression with denominator 




f ^/ V w , (2(ai -a3)(a4-ai)) 

(02 - asjioi - a5J(a4 - 05) W r (02 - 04) + (04 - ai)c2 

(ai -05) ^ 

(5.19) 

Substituting C2 into the expression above we obtain 

2(ai - 02)^(01 - a4)(Q2 - a4)(a3 - 04)^(04 - Q5) 
03-05 

which does not vanish provided that all the a^'s are different. The same 
value of gtt is obtained by taking the limit z — )■ a^ for gu in region p = 0, 
z G (05,04). So we conclude that gtt\p=o is continuous at 2; = 05. A similar 
calculation establishes continuity of gtt\p=o at z = 04; here the relevant 
denominator of the limit z — )■ 04 reads: 

2(02 - 01)^(02 - 04)(a4 - ai)(a4 - 05) . 

The denominator of gu restricted to p = 0, 04 < z < 03 can be written 
as 

2(01 - 02)^(01 - z){a2 - z){z - 05) , 

and is therefore smooth on this interval, extending continuously to the end 
points. 

Non-existence of zeros of the denominator of gu restricted to p = 0, 
03 < 2; < 02 can be proved similarly as for 05 < z < 04. After factorisations 
and cancellations, the numerator of Hx there is a complicated polynomial in 
the aj's, z, and the Cj's, quadratic in C2. The discriminant of this polynomial 
equals 

8(ai - 02)^(01 - 03)^(01 - zf{a2 - 03)^(02 - 05)^(03 - z)(a4 - z){a^ - z) , 

which is negative because of the third-to-last factor. We conclude that gu 
is smooth in a neighbourhood of {p = 0, z G (03,02)}. The continuity of 
gu\p=Q at z = 03 may again be checked by taking left and right limits. 

Non-existence of zeros of the denominator of gu restricted to p = 0, 
a2 < z can again be proved by calculating a discriminant. The numerator 
of Hx there is a quadratic polynomial in C2, with discriminant 

32(ai - 02)^(01 - 03)^(01 - zf{a2 - 04)^(03 - z)(a4 - z){a^ - z) . 

This is negative because each of the three last factors is negative. We con- 
clude that gu is smooth on a neighbourhood of {p = 0, z G (02, 00)}. 
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5.4.3 Ergosurfaces 

The ergosurfaces are defined as the boundaries of the set gu < 0. Their 
intersections with the axis are therefore determined by the set where gu 
vanishes on the axis. We will not undertake a systematic study of those, 
but only make some general comments; see j5] for some results concerning 
this issue. 

Near the points a^ the numerator of gu has the following behaviour: 

~ 4 for ai (see ([2^), 

~ ((02 - a4)(ai - a5)ci + (04 - ai)(a2 - 05)02)^ for 05, 04, 

~ ((02 - a3){ai - 05)01 + (03 - ai)(a2 - 05)02)^ for 03, 

~ ((02 - a3)(ai - 05)01 + (03 - ai)(a2 - 05)02)^(02 - z) near a^, 

~ (2(oi - 03)(o2 - 04) + (oi - 05)0102)^(02 - z) near o^, 

where ~ stands for a manifestly non-vanishing proportionality factor. This 
shows that a component of the ergosurface always intersects the axis at 
z = 02. It also follows from the above that the intersection of the ergosurface 
with the axis {p = 0} contains z = ai and z = 02 as isolated points when 
02 = 0. 

Next, a Mathematica calculation (in which ci has been replaced by its 
values in terms of the Oj's) shows that on (—00, 05) the metric function gu\p=o 
can be written as a rational function with numerator which is quadratic in 
z. Recall that the numerator does not change sign on (—00,05), so gu\p=o 
is continuous with at most two zeros there. But gu\p=o is negative for large 
negative z, while at 2; = 05 we have 

(05 -03) (oi(oi -a5)(o2 - 04) + 02(04 -oi)(o2 - 05)) 
9tt[p = v,z = 05) = 2 

(05 - Oi)(a2 - 05)(05 - 04) (02O1 - O1O2 + 04(02 - Oi)) 

(5.20) 
which is strictly positive. We conclude that gu\p=o always has precisely one 
zero on (—00, 05). 

In Figure 5.1 we show the graph of gu\p=o for a set of simple values of 
parameters. 

5.4.4 gpp and g^^ 

The metric functions gpp = Qzz on p = 0, 2; G (oi, 05) equal 

04 — z 
2(03 - z){z-ari) ' 



(5.21) 
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gtt(p=0, z) 

til- 




Figure 5.1: gtt\p=o as a function of z for oi = 0, 02 = 1, 03 = 3/4, 04 = 1/2, 
05 = 1/4. In this case the ergosurface encloses both horizons. 



and are therefore smooth there. By analyticity, the same expression is vahd 
for z G (—00, as). 

The metric function gpp on p = 0, z £ (05, 04) can be written as a rational 
function of z, with denominator 

4(ai - 02)^(02 - a4)(a2 - z){a3 - 04)^(04 - z){z - as) , 

and is thus smooth near {p = 0,z G (asi04)}|j One checks that for z > 05 
and close to 05 we have 



9pp\p=o 



04 — as 



+ o(i), 



(5.22) 



2(03 -as)(z-as) 
leading to a pole of order one when as is approached from above. Comparing 



with (5.21) one finds that \z — asl x gpp\p=o is continuous at as. 
Next, for z < a4 and close to 04 we have 



9pp\p=o 



as — a4 



2(a3 -a4)(z-a4) 



+ o(i), 



(5.23) 



leading to a pole of order one when 04 is approached from below. 
The metric function gpp on p = 0, z £ (a^, a^) equals 



z-as 



2{z - a3){z - 04) 



(5.24) 



*This denominator has been obtained by substituting the values of k and ci, but not 



C2- 
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with simple poles at 04 and 03. Comparing with (5.23) one finds that 

\z - 04] X gpp\p=o 

is continuous at 04. 

The metric function gpp on p = 0, 2; E (03, 02) can be written as a rational 
function of z, with denominator 

4(ai - 02)^(01 - a3)(a2 - a3)(a3 - 04)^(01 - z){a2 - z){a3 - z){z - 05) , 

which has been obtained by substituting in k, but neither ci nor C2- For 
z > 03 and close to 03 we have 



9pp\p=o 



as -05 



+ o(i), 



(5.25) 



2(03 -a4)(z-a3) 
and there is a first order pole when z = 03 is approached from above. 



Comparing with (5.24) one finds that |z — 03] x gpp\p^Q is continuous at 03. 
Again, for z < a2 and close to 02 we have 

I _ (ai - a3)(Q3 - a5)(2(Q2 - a3)(a2 - 04) + (a2 - a5)cf) 

4(ai - a2)(a2 - asjlos - 04) (02 - z) 

(5.26) 

Since C2 is real, the numerator of the leading term does not vanish. There- 
fore, gpp\p=Q has a first order pole when z = 02 is approached from below. 

The metric function gpp on /) = 0, z G (02,00) can be written as a 
rational function of z, with denominator^ 

4(01 - 02)^(03 - a4)^(-a2 + a4)(a2 - z){as - z){-ar^ + z) . 
Finally, for z > 02 and close to 02 we have 

(ai - a3)(a3 - a5)(2(a2 - a3)(a2 - 04) + (02 - 05)0^) 



gpp\p=o — ,/ w N/ No^ \ + ^(ij • 

4(ai - a2)(a2 - a^jlas - 04)^(02 - z) 

(5.27) 
This coincides with (5.26) except for an overall sign. Again, with C2 being 



real the numerator of the leading term cannot vanish, so the limits from 
above and from below of |z — 02] x gpp\p=o at z = 02 are different from zero, 
and coincide. 
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5.4.5 gtrp and g^^ 

We pass now to the singularities of 

Hy ( Ld^ 

on the axis p = 0. It turns out that the calculations here are very similar 
to those for gu, keeping in mind that the interval (—00,05) was handled in 
Section [2} In particular the lack of zeros of the relevant denominators on 
each subinterval of the z-axis is established in exactly the same way as for 
gu, while continuity at the Oj's is obtained by checking the left and right 
limits. This results most likely from the rewriting 

Fu^ + qFHy 
9tij) — pj^ , 

and noting that, away from the aj's, any infinities of gt^ij\p=Q can only result 
from zeros of FHx- In any case, a Mathematica calculation shows that no 
further infinities in gi^\p=Q arise on the axis from Fu^ + qFHy, and in fact 
the denominators of 5i^|p=o, when this last function is written as a rational 
function of the z's, Oj's, and the Cj's, coincide with those of gtt\p=o- So, we 
find that gtif} is smooth near 

/:={/) = 0, z G (-00, 05) U (05, 04) U (04, as) U (03, 02) U (a2, +00)} . (5.28) 

For the remaining points 02 , . . . , 05 , we write instead 

9t^=gtt{^ + q\ ■ (5.29) 

Using Mathematica we verified that the left and right limits of {oj.^/ Hy)\p=Q 
at aj=i^5^4^3 are equal, but the left and right limit at 02 is not. These are, 
respectively: 



^^^ for ai, 

2(ai— a2)(ai— a4)(a2— 04) 

(02-04) {ai-a5)ci+(a4-ai)(a2-a5)c2 



for as, 04, 



2(ai-a2)(ai-a3)(ffl2-a3) f^^. q- 

(a2-a3)(ai-a5)ci+(a3-ai)(a2-a5)c2 •^' 2 ' 



2(ai-a2)(ai-a3)c2 
2(a-i-az)(a2-a4) + {ai-ar,)c]_C2 ^"^^ "2 



for On 
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(Note that the first hne above contains an inverse power of C2, and so the 



case C2 = requires separate attention; this is handled in Section 5.8.1 ). On 



the other hand, the numerator of gu on p = has aheady been analysed in 



Section 5.4.3, we repeat the formulae for the convenience of the reader 

~c| for ai (see ([2^), 

~ ((02 - a4)(ai - a5)ci + (04 - ai)(a2 - 05)02)^ for as, 04, 

~ ((02 - a3)(ai - 05)ci + (03 - ai)(a2 - 05)02)^ for 03, 

~ ((02 - a3)(ai - 05)01 + (03 - ai)(a2 - a5)c2)^(a2 - 2) near a^, 

~ (2(ai - a3)(a2 - 04) + (ai - 05)0102)^(02 - z) near a^. 

We note that the z-independent terms above all have the same sign when 
C1C2 > 0, hence they are not identically zero. Thus the factors displayed 
here in the numerator of gu can be cancelled with the corresponding factors 



in the denominator in the product gu x {oj^/Hy) arising in (5.29). This 
implies that gt',p\p=Q is continuous for z € M. 
Consider next g^^\p=o-, 

fUJ^ \ Gy 

A Mathematica calculation shows again that the denominator of this func- 
tion, when written as a rational function of z and the a^'s, coincides with 
the denominator of gu\p=o, which has already been shown to have no zeros. 



This, implies that (7^^|p=o is smooth near the set appearing in (5.28). 

From what has been said so far, to prove continuity of g^^ it remains to 
establish continuity of Gy/gu at z = ai. Now, Gy is continuous on p = for 



z gM and vanishes for z > a^ (see Table 5.1) so S',/,^|p=o is continuous at 
{05,04,03,02}. We conclude that g^^ is smooth near the set in (5.28), and 
that g^^\p=o is continuous at all z € M. 

However, the above is not the whole story about g^p^, as we need to know 
where g'^^|p=o vanishes; such points correspond either to lower dimensional 
orbits, or to closed null curves. 

It already follows implicitly from Section [s] that g^^\p=o = for z > 02 
and, in fact, in that interval of z's we have 

<?V^ = 9pp{l + 0{p^))p^ , (5.30) 

as needed for a regular "axis of rotation". This formula is obtained by a 
direct Mathematica calculation, in the spirit of the ones already done in 
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g^il,(P=0, z) 
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Figure 5.2: The graph of g^^ on the axis for ai 
1,^2 = 1. 



0,05 



45^4 



2' "3 



this section. We emphasize that we are not claiming uniformity of the error 
term 0{p'^) above as 02 is approached. 



Note that gpp > away from the axis, and it follows from (5.30) that 
Si/ji/i > for z > a2 and p > small enough. 

The question of the sign of g^^\p=o on the remaining axis intervals is 



addressed in Section [5 . 8 . 3| under the hypothesis that C2 = 0. In Appendix [B] 
we give numerical evidence that g^jp\p=o is positive on {z < 02} for gen- 
eral C2's, see Figure B.2 The values of (7,/„/,|p=o at z = a^ for i = 5,4,3 
can be easily obtained by direct limits computation. As expected from the 
continuity established earlier the right and left limits coincide and are equal 
to 



(a5-a3)(g(ci(ai-a5)(a2-a4)+C2(a4-ai)(a2-a5))+2(ai-a2)(ai~a4)(a2~a4))^ 



(ai -as) (05-02) (05-04) (a2Ci-aiC2+a4(c2-ci))2 

(9(ci(oi-05)(o2-04)+C2(o4-oi)(o2-05))+2(ai-02)(ai-04)(a2-04))^ 
2(01 — 02 )-^ (01— 04) (04— 02) (04— 05) 

(l?(ci(0i-05)(02-03)+C2(03-0i)(02-05))+2(0i-02)(0l-03)(02-03))^ 
2(01 -02)2 (01-03) (03 -a2)(o3 -05) 



for as, 



for 04, 

for 03. 

= a,) > for 



From the ordering of Oj's (A. 15) it follows that g^p^ip = 0,z 
i = 5, 4, 3 if the parameters are distinct. 

Finally, we need to check the signature of the metric. A Mathematica 
calculation shows that near /, as defined in ( 5.28| ), we can write 



detgp, = {f + 0{p^))p^ 



(5.31) 
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where / is an analytic function of z; for example, 



z-°4 7 < n^ ■ 

2{a3-z){z-as) ; ^ ^ "5 , ,^ ^^. 

-^=^ 04 < Z < as . \ ■ > 



2(a3-z){z-a4,) 



(No uniformity near the end points is claimed for the error term in (5.31).) 
The explicit formulae for / on the remaining intervals are too long to be 
usefully cited here. We simply note that we already know that the deter- 
minant of the metric is strictly negative for /) > 0, and thus / < on the 
axis by continuity. However, / could have zeros, which need to be excluded. 



Clearly there are no such zeros in the intervals listed in (5.32). Next, in the 
region z > a2 one finds that / = —h?', where /i is a quadratic function of C2- 
The discriminant of h with respect to C2 reads 

32(ai - 02)^(01 - 03)^(02 - ^4)^(01 - z)'^{a^ - z)(a4 - z){a^ - z) . 

This is strictly negative for z > a2 and we conclude that / does not vanish 
on this interval. 

Taking into account the polar character of the coordinates (p, (/?) and 
(p, ip) near the relevant intervals of z, what has been said so far together 



with formula (5.31) implies that g is a smooth Lorentzian metric on 

M^\{p = 0, zG [05,04] U [03,02]} . 

The missing open intervals, and their end points, need separate attention; 
this will be addressed in Sections 15.51 and 15.61 



5.5 Extensions across Killing horizons 

It is expected that the interval z G [05 , 04] lying on the coordinate axis 
p = 0, corresponds to a ring Killing horizon with topology M x S^ x 5"^, 
while z £ [03,02] corresponds to a spherical Killing horizon, with topology 
M X 5^. The aim of this section is to establish this, modulo possibly the 
end points where the axis meets the Killing horizon; this will be addressed 
in the next section. The construction mimics the corresponding extension 
procedure for the Kerr metric, see also ff] Section 3] or [^. 

Let o G M and let m > be given by 

2 ( ^i ~ '^iX , 2 
m = ^^— — + o , 
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set r± = m i: \^m'^ — a^. As a first step of tfie construction of an extension 
on [ai,aj\ = [05,04] or [ai^aj] = [03,02] we introduce the usual coordinates 
f and 9 for the Kerr metric: 

Ri + Rj ~ -1 I ^j ~ ^i 



+ m, 6l = cos"M^^ -] , (5.33) 

2 \ tti- aj J 

with inverse transformation (see, e.g., [9| (1.133), p. 27]) 

p = y r'^ — 2mf + o^ sin(0) = y (r — r^){f — r+) sin(0) , (5.34) 

z = ^-- — - + {f — m) cos{9) . (5.35) 



Note that in the above conventions we have Oj > o,: 



In the (f, 9) coordinates the flat metric 7 := dp^ + dz'^ remains diagonal, 
\f - mf - [rr? - a^) cos^{9)\ 



(f — m)^ — {m? — a?) 



+ d9^ (5.36) 



^ + (^2 _ ^2) g.^2 q\ ( !i^rf^2 ^ ^^~2 ] (5 37) 

sin^e J \ P^ j 

= R^Rj (jz %z r + d~9A (5.38) 

= R.R, ( ^ (^-^-- + ^l' "-^^ V '^ - ^') df^ + A , (5.39) 
y p [o-i aj) J 

where the various forms of the metric 7 have been listed for future reference. 

The essential parameter above is m^ — o^, in the sense that a change of 
m and o that keeps m^ — o^ fixed can be compensated by a translation in f, 
without changing the explicit form of 7. The replacement of ^/Trfi — a?' by 
—^/•w? — a^ can be compensated by a change of the sign of (f — m), which 
again does not change the explicit form of 7. 

We have, near p = 0, for Ui < z < aj, with error terms not necessarily 
uniform over compact sets of z, 

= 4(a. - .)V, - .)% ) _ p,,„) 

/)^(Oj - Oj)^ 

7,-,- = \{z-ai){z-aj)\+Oip^). (5.41) 

Now, the Black Saturn metric depends upon p through p^ only, with the 
latter being an analytic function of f and 9. In the new coordinate system all 
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the metric functions extend analytically across {p = 0,z G {ai,aj)} except 
gfr, which has a first order pole in f at f = r-t. In the original coordinate 
system we start with f > r_|_ and it is not clear whether or not r = r- can 
be reached in the analytic extension, but we need to get rid of the pole at 
f = r_|_ in any case. For this, it is convenient to continue with a general 
discussion. We consider a coordinate system {x'^,y) = (a;'^,x*) = {x^,x^,y), 
where p runs from to n — 1, and we suppose that: 

1. The metric functions g^i, are defined and real analytic near y = yo, 
except for gyy which is meromorphic with a pole of order one at yo- 

2. The determinant of the metric is bounded away from zero near y = yo- 

3. There exists a Killing vector field ^ of the form 

C = do + a'di , 
for some set of constants a*, such that all the functions 

vanish at y = yo. 
In our case the first condition has just been verified with 

y = r , yo = r± . 

The determinant condition holds by inspection of the metric, see Ta- 
bles O and O 

The third condition is verified by a Mathematica calculation, leading 
to a Killing vector dt + 0,^3 d^ , where 

^ __f 2(ai - a2)(ai - Q3)(Q2 - Q3) V 

^ V("2 -a3)(ai -a5)ci + (as -ai)(a2 -a5)c2 / 

satisfying the condition on (03,02), and the Killing vector dt + ^^51x529^ , 
with 

^ __f 2(ai - a-2){ai - ai){a2 - Q4) V 

^ ^^ V("2 -a4)(ai -a5)ci + (a4 -ai)(a2 -a5)c2 / 

satisfying the condition on (05,04). A rather lengthy MATHEMATICA calcu- 
lation shows that the O's are finite for distinct a,'s. 
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We introduce new coordinates (x^,y) = {x'^,x ,y) = (x'^,x*) by the 
formula 

x° = x\ x' = x'- a'x" . (5.42) 

This coordinate transformation has Jacobian one. Writing gf^i, for g{dxi^ , S^m ) , 
our hypotheses imply that we can write 

% = (y - yo)xf. , 9yy = (^3^ ' (^-^^^ 

for some functions Xfn h, all analytic near i/q. 

Since the metric functions are now independent of x'^, the next coordinate 
transformation 

dx° = dx^ + f{y)dy , 5^ = x^ , y = y , 

again with Jacobian one, does not affect the analyticity properties of the 
functions involved. We have 

Sooidx^f + gyydy'^ = {y - yo)Xd{dx° - fdy) + -j—^ — zdy^ 

= (.y-yo)Xo(.dx'^f -2{y-yo)xofdx°dy 
^ h+{y-y,fx,f 

[y - yo) 

AC := — lim — 

y^yo Xq 

is a positive constant. Keeping in mind that Xq is negative while h is positive, 
and choosing / as 



Assume that 



/ = ^^ , (5.45) 

y-yo 

one obtains a smooth analytic extension of the metric through y = yo, since 



then the singularity in (5.44) is removable; similarly 

g^-dx^dx' = {y-yQ)Xi{d5p - fdy)dx' 

= {y- yo)Xidx^dx' - X'iVi^dydx' . 

The determinant of the metric in the coordinate system x^ equals that in the 
original coordinates, and so the extended metric is Lorentzian near y = yo- 

It remains to show that this procedure applies to the BS metric, with 
x^ =t , y - yo := r - r+ , (x^) = {ip, ip, 6) , 
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where f and 6 have been defined in (5.33). We have 



(aj — aAp^ ^ , A. 

A{aj - z){z -Ui) 

hence 

(f - r_) sin^ (9 = -^-^. ^^^- — ^ + 0{p^) , 

(Oj - ai) 

with the error term not uniform in z near the end points. On (05, 04) or on 
(03,02) one needs to calculate the limits 

h\r=r+ = lim— "^ . ~ X Iini(/3^7f=f=) . 

p~!>o (r - r_)sm''6' p^o 

Letting J] = ^151x52 on (05, 04), respectively 0, = Clga on (03, 02), one further 
needs 

Xo\r=r+ = lim {p^'^g{dt + 0,8^, dt + i^5^)(r - r_) sin^ ^i) . 

A surprisingly involved Mathematica calculation shows that at p = the 
quotient h/xo equals, up to sign, 

(04 - a5)(2(Qi - a2)(ai - Q4)(Q2 - Q4) + ((a2 - Q4)(qi - Q5)ci + (04 - ai)(a2 - Q5)c2)g)^ 

8(ai - 02)^(02 - a4)(a3 - 04)^(04 - ai) 

on (05, 04), and 

(as - Q5)(2(ai - a2)(ai - Q3)(Q2 - Q3) + ((^2 - Q3)(qi - 0-5)01 + (03 - Qi)(a2 - Q5)c2)g)^ 

8(ai - 02)^(02 - a3)(a3 - ai)(a3 - 04)^ 

on (03,02). As those limits are constants, we have verified that, within the 
current range of parameters, the Black Saturn metric can be extended across 
two non-degenerate Killing horizons. 

5.6 Intersections of axes of rotations and horizons 



It follows from (5.33) that 



Ri = r-r+ + "^ °^(cos6i+l), (5.46) 

Rj = f-r+ + "^'~"' (l-cosg) , (5.47) 

Hi = (f-r+)(l-cos6i) , (5.48) 

Hj = (f-r_)(l-cos6i) , (5.49) 
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so that //j, ^j, Ri and Rj are smooth functions off and cos^rl Furthermore, 



it follows from (5.34) that the function p^ is a smooth func tion of f and of 
slv? 9 = 1 — cos'^ 6, similarly z is smoo t h in c os 6 by (5.35[), which implies 



that the remaining /x^'s (compare ( 5.51 )-( 5.52)) are smooth in f and cosO. 



Now, consider any rational function, say W , of the /ij's and p^, which is 
bounded near f = r_|_, = 0. Boundedness implies that any overall factors 
of f — r+ in the denominator of W are cancelled out by a corresponding 
overall factor in the numerator, leaving behind a denominator d(f, 6) which 
can be written in the form 

d{f, 6) = /(cos 9) + {r — r-^)g{r, cos 9) , 

for some functions / and g which are smooth in their respective arguments. 

// 

d(f = r+,0) = /(l) 

does not vanish at 9 = 0, then the denominator d is bounded away from 
zero near r = rj^ and 9 = 0. This in turn implies that 1/d is smooth in a 
neighbourhood of the point concerned, and therefore so is W. 

An identical argument applies at 9 = tt. 

This reasoning does not seem to apply to uj^ , because of the square roots 



there. However, as mentioned in Appendix A.l these appear in the form 



MqMi M0M2 /M1M4 /M2M4 



X V ^x V '-^x V ^x 

One checks that the expressions under the square root are squares of rational 
functions of the //^'s, and of p'^, and so the metric functions involving lo^ are 
also rational functions of the /ij's and p^. 

Since we have already shown that the suitably reduced denominators of 
all the scalar products g{X, Y), where X,Y £ {dt, d^, d^^}, have no zeros at 
the axis points p = 0, z = ai, we conclude that the corresponding metric 
coefficients are analytically extendible, by allowing f to become smaller than 
r+, including near the intersections of axes of rotation with the Killing 
horizons. 

One similarly establishes analytic extendibility of g^^j : 

_ {gtt + 2gtv.^ + ffv,^r?^)VK 



It should be kept in mind that cos^ is a smooth function on the sphere, but sin 6 is 
not. 
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Here we have already verified that gtt + 2gtip^+ g^^^'^ is an analytic function 
of f and cos^, and extendibility of g^i readily follows from the fact that Vt 
has been chosen so that this function vanishes at f = rj^. 
Finally, gyy is given by the formula 

\/^%-+(^-^+)5r>- 
9yy = f_r+ ' ^^'^"^ 

To analyse this metric function, by a Mathematica calculation we verified 
that the reduced denominator of (f — r^)gfr does not vanish at f = r_|_, 
and hence this function extends across f = r^ as an analytic function of f 
and cos^. Keeping in mind that the same has already been established for 



y/Rgiy, we find that the numerator of (5.50) extends across f = r+ as an 
analytic function of f and cos 9. Analytic extendibility of gyy follows again 
from standard factorisation properties of such functions. 
We next analyse g^^ near p = 0, z = a^. Now, 

and we need to understand the behaviour of the functions above near r = r^, 
6 G {0, vr}. For i ^ 5 we have 

fJ-efJ-b + P^ = Ur -r^)sm^ 9 + iJ.e{l- cos 6)) {f -r+) , (5.51) 
and since 

2 2 

/il = - — ^ —-^ r (5.52) 

Ri + z — ai 2(04 — ai) 

near /O = 0, z = 04, for £ = 2, 3 we can write 
1 ( ~ i^l 



1 + cos / 



1 I ~ l^e 
fJ-efi-b + P = \r -r- -\ 

V 1 + cos / 

2 2(f — r_) N . 2 /I 
P4P5 + P = jr{r -r+)sm 9 , 

1 + cos 9 

, 2 _ f "■^~ ^^ i -i\ 2 
piUe + p w Kl p 

V 04 - ai J 



{f — r^) sin 9 , 


(5.53) 


(f — r+) sin^6' , 


(5.54) 


9, 


(5.55) 



04 — ai 



-(f — r_)(f — r_|_) sin 9. (5.56) 



Finally, for £ = 1,4, 5, 

p,iPi + p ss p = (f — r_)(f — r+) sin 9. (5.57) 
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Encoding this behaviour into a Mathematica calculation, one finds that 
Qgg is uniformly bounded in a neighbourhood of r = r_|_, cos 9 £ {±1}, with 
non-vanishing value of the denominator as needed above. This establishes 
smoothness. Similarly g^^/ slv?6 is smooth near those points. 

Now, away from, and near to, the event horizons, the map (p, z) i— )■ 
(f, 9) is a smooth coordinate transformation. From what has been already 
established, the two-dimensional metric 

g^gdO'^ + g^^dip^ (5.58) 

is thus a smooth metric for f > rj^, f close to r_|_, in particular there is 
no conical singularity at the rotation axis for d^p in this region. But the 
arguments just given show that this metric extends smoothly across f = r_|_, 
which finishes the proof of smoothness of the whole metric up-to-and-beyond 
the horizon near f = r-|_, ^ = 0. 

A similar analysis applies near 05, 03 and 02; in this last case, one 
considers the two-dimensional metric 



g^Qde"^ + g^ij^d^p^ 



instead of (5.58) 



5.7 Event horizons 

Consider the manifold, say ^ , obtained by adding to the region f > r_|_ those 
points in the region r^ < r for which the metric is smooth and Lorentzian. 
Then the region r_ < f < r-(_ is contained in a black hole region in the 
extended space-time, which can be seen as follows: Note, first, that 5^^ 
vanishes at J^ := {f = r_|_} = {y = yo}, which shows that J^ is the 
union of two null hyper surf aces. On each connected component of J^ the 
corresponding Killing vector X = dt + ^d^ is timelike future pointing for 
y > yo close to yo, and so by continuity X is future pointing on Jif. This 
implies that J^ is locally achronal in the extended space-time: if a future 
directed timelike curves crosses Jif through a point p G J^, it does so 
towards that side of TpJ^ which contains the component of the set of causal 
vectors at p containing X. Since J^ is a (closed) separating hypersurface in 
^, this implies that any timelike curve can cross J^ only once. From what 
has been said it follows that the region r_ < f < r+ is contained in a black 
hole region of {^,g). 

In particular we have shown that the black hole region is not empty. A 
standard argument (compare [!2,, Section 4.1]) shows that J^ coincides with 
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the black hole event horizon in ^. Note that this is true independently of 
stable causality of (^, g), or of stable causality of the d.o.c. in (^, g). 

Some more work is required to add the bifurcation surface of the horizon, 
a general procedure how to do this is described in [11 J. 

5.8 The analysis for C2 = 

We turn our attention now to the Black Saturn solutions with C2 = 0, where 
the formulae simplify sufficiently to allow a proof of stable causality of the 
d.o.c. 



First note that (4.2) implies that the condition C2 = leads to ci 7^ as 
the only restriction on ci. However, it implies a fine-tuning of the parameters 
ai. One may easily check that the minus sign solution for C2 cannot vanish 



if the ordering (A. 15) of the Oj's is assumed. However the plus sign solution 
may lead to the vanishing C2 under certain additional conditions. Namely 
the resulting equation 

\/(a3 - ai)(o2 - a4)(a2 - a5)(a3 - «5) = (0-2 - «i)(a3 - 04) , 

quadratic in 05, may always be solved for 05 = 05(01,02,03504) £ ^', the 
condition that < 05 < 04 is then equivalent to 

04 < (02 + 0103 - 2o2a3)/(ai - 03) . (5.59) 



This is more transparent in terms of the variables Ki G [0, 1] defined by (5.3), 



as then (5.59) becomes 

/€! > -^— , (5.60) 

2 - K2 



see Figure 5.3 In the further analysis one should keep in mind that 05 is no 
more an independent parameter. 

Notice that C2 = implies q = and k = 1. 

5.8.1 Smoothness at the axis 



Smoothness of the Black Saturn solution for p > 0, proved in Section 5.3 
holds also for the C2 = case, hence only the analysis on the axis of rotation 
needs separate attention. We shall proceed in the same way as in Section 
53 

We start with an analysis of the behaviour oi g^^ on the axis. For z < ai 
it may be written as a rational function 

2(ai - 03)^(02 - z){z - a2){z - ai){z - 05) + ci^(ai - 02)^(01 - 05)^(03 - z 
(oi - 03)2(01 - z){z - a2){z - 04) 



5 THE ANALYSIS 



34 




«i 



Figure 5.3: The variable ni runs along the horizontal axis, while K2 runs 



along the vertical one. The inequality (5.60) corresponds to the shaded 
region. 



To avoid the singularity at z = ai we need to fix ci as to have a finite limit. 
Miraculously, this condition leads to the same formula ci as obtained in 
section 2 for C2 7^ 0. This is somewhat unexpected, since we have set C2 to 
zero as an alternative to fixing ci. With this choice of ci regularity on the 
axis of many metric functions has already been established, and we would 
be done if not for the fact that some of the formulae derived so far involve 
explicit inverse powers of C2. So it is necessary to repeat the analysis at the 
axis from scratch. 

Several formulae are much simpler now. For instance, one checks that 
in the region oi < z < 05 on the axis g^^ is given by the same formula 
as for z < ai. Hence we conclude, that g^^ is smooth and bounded for 
{p = 0, z < as}. 

In the subsequent axis interval, 05 < z < 04, g^^ is a rational function 
with denominator 



2(ai - 02)^(01 - a4)^(a3 - z){a5 - z) - ci^(ai - 05)^(02 - zf{ai - z) , 

which cannot vanish, being a sum of two negative terms. At both end 
points of the investigated interval one of the terms in non-zero, which shows 
boundedness. 
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Moving further to the right we obtain a simple formula for 5^^: 

2(ai - z){a2 - z) 



(5.61) 



(as - z) 

which immediately implies continuity for 04 < z < 03. We note that this is 
strictly positive, and therefore near that axis interval g^^ is strictly positive 
as well. 

In the region 03 < z < 02 the denominator of g^^ is more complicated: 

(ai - 05)^01(02 - zf{a:>, - z) - 2(ai - 02)^(01 - 03)^(04 - z){a^ - z), 

but does not vanish, being a strictly negative sum of two non-positive terms. 

In the region z > 02 for vanishing C2 the function g^^ is proportional to 
q'^. Since C2 = implies g = 0, we conclude that g^^ vanishes for z > 02, as 
already seen for general values of C2 in any case. 

The analysis of gu is similar. For p = and z < 05 the metric function 
gtt is a simple rational function, 

(«i - z)iz-a3) 
[z - a2){z - a4) 

which is clearly continuous in the region z < a^. For 05 < z < 04 the 
denominator of gu reads 

(ai - 05)^0^(02 - z)'^{a4 - z) + 2{ai - 02)^(01 - 04)^(03 - z){z- 05) . (5.63) 

with both terms manifestly positive in the region 05 < 2; < 04. We conclude 
that gu is smooth on 05 < z < 04, bounded on 05 < z < 04. 
Next, for 04 < z < 03 the denominator of gu reads 

2(ai - 02)^(01 - z){z - a2)iz - 05) , 

thus it cannot vanish for 04 < z < 03. Moving further to the right we find 
the denominator of gu 

(ai - a5)^Ci(a2 - zf{a3 - z) + 2(ai - 02)^(01 - 03)^(04 - z){z - 05) 

as a sum of manifestly negative terms on 03 < z < 02. Also the end points 
are singularity-free. Finally, for z > 02 gu equals 

i^""ii^""i , (5.64) 

(ai - z)[z - 03) 

hence it is continuous. This proves directly absence of singularities for gu 
on the axis in the case of vanishing C2- 

The analysis of gf^ can be carried out along the same lines, and is omit- 
ted. 
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5.8.2 Causality away from the axis 

We have not been able to estabhsh non-existence of closed timelike curves for 
a general Black Saturn solution, though we failed to find any in a numerical 
search, see Appendix [B] However, if one imposes the condition C2 = the 
metric formulas simplify sufficiently to allow a direct analysis. Indeed, the 
explicit formula for g^^ in the case of vanishing C2 (and consequently q = 0) 
reads 

p2 (ci^Mi + Mo) (Mo//i2 - ci2Mip2) "• sRI ■ 

Outside the axis (p > 0) the ordering of //j's is the same as those of aj's 
and all the functions Mi are strictly positive. Both the numerator and 
denominator of gf^^ can be regarded as quadratic functions of cf. Let us 
first investigate the possible zeros of the denominator: 

5(ci) - ^ ci - ^^ or ci - ^^^2 . 

Clearly only the second one is relevant since the first one would lead to an 
imaginary coefficient ci. On the other hand, the equation /(cf ) = has two 
solutions: 

^ Mo {2Ri (Ri ± V^i' - P^) - P' 



c± 



Mi/)2 

To make this result more transparent let us express Ri in terms of /Ui and p 
Then c^ may be written as 



c± 



Mip2 2pl 



From the explicit form of /zi one can easily see that sign(/uf — p2) = sign(ai - 
z), thus we have: 

2 ^Wo/^ 2 ^OP' 

tor z < ai c_ = — - — ^ , c, — 



Mip2 ' "+ Ml/if ' 

for z > ai c2 - ^^ ^2 _ M)pi 
tor z > ai c_ - ^^^2 , c+ - ^^^^ • 
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We see that in both regions one of the zeros c^ of the numerator cancels 
the zero of the denominator, which provides an alternative explicit proof of 
regularity of g^^ for p > 0. Moreover we find: 



S'V"/' ~ „2 ('„ 2 )i,r > i\/r \ \ HJ ..2 '^1 



Keeping in mind that the parameter ci has been fixed to guarantee the 
regularity on the axis, to obtain a sign for g^^ for p > it remains to show 
that the equality 



ci 



Mop2 



Ml/if 
can never be satisfied away from the axis. For this, we shall make use of the 



formula (5.8 ) expressing cf in terms of /ij's. By subtracting the two formulae 
for cf we obtain 



[l-n^if^l - /"2)^(/Ul - /i4)(/"l - P5) (/"l/"3 + P^) 

+ (/il - ;U3) (/il;U2 + p^) (/il/U4 + p^) (/il/U5 + p2) j = . 

The overall multiplicative coefficient in the first line is strictly negative, 
whereas the term in parenthesis across the second and third lines is a poly- 
nomial in p with coefficients that can be written in the following, manifestly 
negative form 

/^l ((/^l - At2)^/U3(/^l - /^4) + P5 (/^i(/U4 - Ps) + (pi - /^4)^3(^2 + (/i2 - Ml)))) 
+p2pf (^P2(m - /U3)(2/i4/U5 + M2(At5 + f^l)) 

+ (P4 - Ml) (m2(M5 - Ms) + Mi((M1 - M2) - M2)(M5 - Ml)) ) 

+p'^Mi(mi - Ms) (m2 + M4M5 + 2m2(M4 + Ms)) 
+p'^Mi(Mi - Ms)(2m2 + M4 + Ms) 
+P®(Mi - Ms) • 

It follows that g^^ > for p > when C2 = 0. 

It turns out that an alternative simpler argument for positivity can be 



given as follows: Using (5.8) we may write g^^ in terms of pi and p. The 
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functions ^j satisfy the same ordering as a, (A. 15) (see (5.2)). The strict 



version of the ordering (A. 15) impUes a strict ordering of the ^j's for p > 0. 
Assuming that, we may make the positivity of g^^ exphcit by expressing it 
in terms of the positive functions 

A51 = /U5 - /^i, A45 = /X4 - Ai5, A34 = /U3 - //4, and A23 = /i2 - fJ'3- (5.65) 

The numerator and denominator of g^^ are polynomials in Aij, /ii and 
p, the explicit form of which is too long to be usefully exhibited here. By 
inspection one finds that all coefficient of these polynomials are positive, and 
since the Ajj's, pi and p are positive, both the numerator and denominator 
of g^^ are positive. 

5.8.3 Causality on the axis 

We turn now our attention to the axis. By continuity, we know that g^^ 
at p = is non- negative. It therefore suffices to exclude zeros of 5'^^|p=o- 
Equivalently, whenever we find a manifestly non-zero value of ^^^(0, z), we 
know that this value cannot be negative. 

Now, at p = and for 2; < ai we replace z hy w := z — ai < 0, and find 
that g^^ there is a rational function with denominator 

(ai - a3)(ai - 02 + w){ai - a^ + w) , 

which is seen to be strictly negative for w < 0. On the other hand, the 
numerator is a third-order polynomial in w: 

2{a2 - ai) X (Sal - 01(02 + 2(203 + 04 + 05)) 

+ 01(20203 + 303(04 + 05) + 0405) — 02(03(04 + 05) — 0405) — 2030405J 

+2zi;(o3 - oi) (6oi(oi - 02) - 301(04 + 05) + 02 + 202(04 + 05) + 04O5) 
+2w (03 — oi)(4oi — 2o2 — 04 — 05) 
+2u;^(o3 -01) . 

Unless explicitly indicated otherwise, the remaining analysis uses the choice 
of origin and scale given by oi = and 02 = 1, which involves no loss of 
generality for checking the sign of g^^ . The above reduces then to 

-2((2o3 - 1)0405 + 03(04 + 05)) + 2o3u;(o405 + 2(04 + 05) + 1) 
— 203U; (04 + 05 + 2) + 203U; . 
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Each monomial in the above polynomial is manifestly strictly negative for 
w < 0, except perhaps for the zero-order term. However, when C2 = 0, in the 



current choice of scale we necessarily have 03 > 1/2 by (5.59), which makes 
manifest the negativity of the zero-order term as well. Hence fl'^^|p=o > 
for z < ai. 

The interval (ai, a^) requires more work, and will be analysed at the end 
of this section. 

For z £ (05,04) we obtain 

204(2 — l)(a3 — z) 
''~° 03 {a^ia^iz - 2) + 1) - 05(2 - 1)2) + 04(05 - z) ' 

which has no zeros in [05, 04], and thus is positive there. 



Positivity on [04,03] follows already from (5.61). 
For z G (03,02) we obtain 

_ 203(2;- l)(a4 -z) 



03(04052: — 204O5 + 04 + 05 — 2) — 0405(2 — 1) 



which again has no zeros in [03, 02], and hence is positive there. 

We already know that {p = 0,z > 02} is a regular axis of rotation for 
d^, so there are no causality violations there associated with d^. 

We consider now the interval (01,05) = (0,05). There we find 

I f 

gii)ip\p=i} 



03(-l + 2)(-a4 + 2) 

with 

/ := (03 (04 (-(05 + 2)2 + 205 + 2^ + 1) + (2 - 1)^(05 - 2)) - 0405) . 

Suppose that there exists 2 in this interval such that / vanishes for some 
< 05 < 04 < 03 < 1. Since / does not change sign, this can only occur if 
at this value of 2 we also have 



daj = daj = daj = . 



Now, 



daJ = 2(-a5 + 03(-a5(-2 + 2) + (-1 + 2)^)) , 
daJ = 2(-a4 + 03(-a4(-2 + 2) + (-1 + zf)) . 
The resultant of these two polynomials in z is 

16(03 - l)^a|(o4 - 05)^ , 
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which is strictly positive in the region of interest, hence g^^ is also strictly 
positive on {p = 0, z £ (01,05)}. 

An alternative argument for positivity at p = can be given as follows: 
Since all terms in the numerator and denominator are non-negative one 
needs to check zeros of the numerator and denominator. The analysis is 
done separately on each interval (cj, Uj). Before passing to the limit p = 0, 



for z > Ui the functions Ajj (as defined in (5.65), and which necessarily 
vanish at p = 0) are replaced by positive functions Ajj such that Ajj = 
p'^Aij. Furthermore we introduce pi = p^pi for z > ai. Substituting 
these expressions in respective intervals of z, cancelling common factors and 
taking the limit p — t- one obtains expressions for the numerator and the 
denominator of g^^ at p = 0. These expressions turn out to be polynomials 
with all coefficients positive. For example for z G (134,133) we obtain the 
manifestly positive expressions 

I _ (A23 + A34)(A51 + Al)(l + (A23 + A34)Al)^ 

^'^'^''=°" /ii(l + (A23 + A34)Al)^ 

and for z G (03, 02) 

A23(A34 + A45 + A5i)(A5i + pi)il + A23Al)^ 



9ipil)\p=o 



/il(A45 + A51 + A34(l + A23/il)^ + A23(A45 + A5i)/ii(2 + A23(A5i + pi))) 



It turns out that the denominator never vanishes and the numerator van- 
ishes, as expected, only at the axis of rotation of d^ (z > 02). 

5.8.4 Stable causality 



Using (A.l) 



g{Vt,Vt)=g' 



Gy 



we conclude from what has been said so far and from Table 15.11 that t is a 
time-function on 

{p>0}U{p = 0,z^ [as, 04] U [03, 02]} , (5.66) 

except perhaps for p = 0, z > a2. There we find 

p-s>o p"^ 2{z - a2){z - a^) 



which ends the proof of stable causality of the region (5.66) when 02 = 0. 



(The blow-up at z = 02 appears surprising at first sight, but turns out to 
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be compatible with a smooth axis of rotation, as clarified in Section 5.6 



compare also (5.30).) 



A The metric 

A.l The metric coefRcients 

The Black Saturn line element |4| reads 

H,, r .. /Wi/j 



ds 



dt+r-f+q]d^ 
tlx L ^rLni 



+ HA k'p(dp^ + dz^) +^di^^ + ^ V 



Hy 



Hx 



where /c, q are real constants. The contravariant components of the metric 
tensor are g^^ = Hy/iH^Gy), gPP = g'' = l/g^p, g^^ = 1/g^^ and 



tt _ H^ Hy fu}^ 



-^y -^x^y \-^y 



JhMj_ tip ^ Hy f^^P^ 

&y HxCry \Hy 



If we let 



Ati := \/p2 + {z - aiY - {z- Oj) , 
where the aj's are real constants, then 



G, 






(A.l) 



(A.2) 



P = {fJ'3f^4 + P ) il^l /^5 + /O )(A*4 ^5 + /O 



(A.3) 



Hx = F 



-1 



Mo + cl Ml + clM2 + ci C2 M3 + clcl Mi 



H„ = F 



-1 /f3 



M,^-clM^^- 

P2 fJ'l fJ'2 



(A.4) 



C2 M2 ^ + Cl C2 M3 + C1C2 M4 — 

r Pi 



(A.5) 
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where ci and C2 are real constants, and 

Mq = ;U2Ai5(^l -/X3)^(/U2 -/^4)^(/0^ + Ml^2)^(p^ +Ail/U4)^(/0^ +/U2/U3)^, 

(A.6) 

Ml = /uf /i2/^3/W4/^5/0^(/Wl -^2)^(/i2 -/^4)^(m -/U5)^(P^ +/^2At3)^, 

(A.7) 

M2 = M2/^3/«4M5P^(/"l -/i2)^(/«l -/i3)^(P^ + W/U4)^(P^ +/"2M5)^, 

(A.8) 

Ms = 2/ii//2^3/"4M5(/"l - /U3)(/"l -/«5)(^2 -/"4)(/0^ +^?)(P^ +^i) 

X (p2 + n^ ;^4)(p2 ^ ^2 /i3)(p^ + /U2 /is) , (A.9) 

M4 = l4l^2l4l4ilJ-l- l^^fiP^ + l^ll^2fiP^ + 1^2 1^5^ , (A. 10) 

and 

^ = /il/"5(m -/"3)^(/^2 -/^4)^(p^ + /ilAt3)(P^ + /^2/i3)(P^ + /il/i4) 

5 

X (p^ + /i2 /U4)(/5^ + /i2 /U5)(P^ + /i3 Ms) JJip^ + /i' ) • (A.ll) 

1=1 

Furthermore, 

G, = ^^, (A.12) 

/i4 

and the off-diagonal part of the metric is governed by 

ci Ri VMoMi - C2 R2 VM0M2 + c? C2 R2 VM1M4 - ci ci Ri VM2M4 



w^ 



(A.13) 



Here i?i = y^p^ + (z — aj)^. We note that the square roots in (A.13) are an 
artifact, in the sense that the functions 

MqMi M0M2 M1M4 M2M4 

^x ^x ^x ^x 

can be checked to be complete squares, which implies that their square roots 
can be rewritten as rational functions of the jUj's, p^, and of the free constants 
appearing in the metric. 

The determinant of the metric reads 

detg^^ = -p'Hlk^P^ . (A.14) 
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A. 2 The parameters 

Here we summarise the restrictions imposed in [41 on various parameters 
appearing in the metric. The parameters Oj are ordered as 

fli < ^5 < ^4 < «3 < ^2 , (A. 15) 

but throughout this paper we assume that the inequahties are strict. 
Boundedness of gu near ai leads either to C2 = or to 



ci 



'2(a3 -ai)(a4 - ai) 
as - ai 



(A.16) 



This last condition follows also from the requirement of boundedness of g^^ 
near ai when C2 = 0, and thus (A.16) needs to be imposed in all cases. A 



choice of orientation of ip leads to the plus sign. 



From Table 5.2 continuity of the metric at {p = 0, z < ai} leads to the 
condition 

^ ^ 2(ai -a3)(Q2-a4) .^ -^^. 

2(ai - a3)(a2 - 04) + (ai - «5)ciC2 ' 

which can be checked to be finite when the value of C1C2 is inserted. 
Asymptotic flatness requires 

2C2K1 

q = , 

2ki - 2kiK2 + C1C2K3 



as well as 



where 



2«;i(-l + K2) 



Y^(-2ki(-1 + K2) + CiC2K3)2 

Qi+2 — Ql 

Kj : — , 

a2 — ai 



which can be checked to be consistent with (A. 17). 

A conical singularity on the rotation axes of d^p is avoided if 



C2 



\/2( 



04—02 



±(qi - Q2)(Q3 - Q4) + \/(ai - Q3)(Q4 - Q2)(«2 - Q5)(Q3 - Q5) 
i/(ai - a4)(a2 - a4)(ai - a5)(a2 - a5)(a3 - 05) 
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B Numerical evidence for stable causality 

In this Appendix we present numerical results that support the conjecture 
that g^p^ is positive away from points where d^ vanishes. Regions where 
g^^ vanishes or becomes negative contain closed causal curves. On the 
other hand, the conjecture implies stable causality of the domain of outer 



communications, see Section 5.8.4 



While our numerical analysis indicates very strongly that g^^ is never 
negative in the region of parameters of interest, it should be recognized that 
the evidence that we provide concerning null orbits of d.^ is less compelling. 

The metric component g^^ is a complicated function of p, z and the five 
parameters aj=i^,,.^o- This function is sufficiently complicated in the general 
case that there appears to be little hope to prove non-negativity analytically. 
We gave a complete analytic solution of the problem in Section |5.8| only for 
C2 = 0. In general, we turn to numerical analysis. The idea is to find an 
absolute minimum of g^^ . 

The original phase-space of this minimization problem is seven dimen- 
sional. One may use translation symmetry of Black Saturn solution to re- 
duce the dimension by one. We do this via the choice oi = 0. Next choosing 
05 — oi as a length unit leads us to a five dimensional minimization problem. 
Our five variables are p, z, ^45, ^34, ^23, where dij = ai — aj. All of them 
are real and in addition p > 0, dij > 0. 

The minimization procedure starts at a random initial point and goes 
towards smaller values of g^^. For general p > we use an algorithm with 
gradient — the so called Fletcher-Reeves conjugate gradient algorithm. The 
limit /) — >• is non-trivial, therefore it has to be studied separately. In this 
case, the values of the metric functions are given by different formulas for 
different ranges of z coordinate. The expressions for the gradients are huge 
and we did not succeed in compiling a C++ code with these definitions. 
Therefore, for p = we use the Simplex algorithm of Nelder and Mead. 
This algorithm does not require gradients. Both algorithms are provided by 
the GNU Scientific Library [6j. 

The minimisation procedure stops when the computer has attained a 
local minimum by comparing with values at nearby points, or when the 
minimizing sequence of points reaches the boundary of the minimization re- 
gion (coalescing Oj's). All local minima found by the computer were located 
very near the axis p = 0, where the results were unreliable because of the 
numerical errors arising from the divisions of two very small numbers, and it 
is tempting to conjecture that g^^ has non-vanishing gradient with respect 
to {p, z, tti) away from the axis, but we have not able to prove that. 
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The numerical artefacts, just described, were filtered out as follows: Each 
value of g^^ at a local minimum, as claimed by the C++ minimisation pro- 
cedure was recalculated in Mathematica. If the relative error was bigger 
than 10~^, then the point was classified as unreliable and excluded from the 
data. In particular all points at which C++ claimed a negative value of g.^^ 
were found to be unreliable according to this criterion. 



Figure B.l illustrates a roughly quadratic lower bound on 



9i)ip\p>o,ze[-z: 



max i^max 



] ' 



with a slope depending on the collection (zmax , dij] 




Figure B.l: The values of g^^ as a function of p at the end of the minimiza- 
tion procedure; this occurs either at local minima, or at points where the 
minimizing sequence leads to coalescing Oj's. The three samples a), b), c) are 
presented with different grey intensity (from low to high, respectively). The 
initial parameters (z, dij) for the minimization procedure were randomly cho- 
sen, uniformly distributed in the intervals a) z £ (—150,301), dij £ (0,50), 
b) z G (-150,226), d^j e (0,25), c) z e (-150,166), dij G (0,5). For each 
sample, the minimum of g^^ is proportional to p^. 



In Figure 



B.2 



one observes a linear lower bound on 5^,/,|p=o for z < ai, 
with a slope approximatively equal to —2 with our choice of scale 05-01 = 1. 
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-150 -100 



Figure B.2: The values of g^^ for p = at the end of the minimization pro- 
cedure; this occurs at points where the minimizing sequence leads to coalesc- 
ing Oi's. The initial parameters {z,dij) for the minimization procedure were 
randomly chosen, uniformly distributed in the intervals z G ( — 150,301), 
d,j G (0,50). 



The numerical results presented in this section support the hypothesis 
that g^^ is never negative in the region of parameters of interest, vanishing 
only on the axis of rotation {p = , z > a2}- 
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